Spin domains in ground state spinor Bose-Einstein condensates 
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^ ' Bose Einstein condensates of dilute atomic gases, 
I— characterized by a macroscopic population of the 
Q^juantum mechanical ground state, are a new, weakly 
interacting quantum fluid jlj |^]. In most exper- 
lents condensates in a single weak field seeking 
^tate are magnetically trapped. These condensates 
ddan be described by a scalar order parameter similar 
£So the spinless superfluid 4 He. Even though alkali 
j la toms have angular momentum, the spin orientation 
0j£ not a degree of freedom because spin flips lead 
"^o untrapped states and are therefore a loss pro- 
cess. In contrast, the recently realized optical trap 
^gor sodium condensates Q confines atoms indepen- 
dently of their spin orientation. This opens the pos- 
sibility to study spinor condensates which represent 
system with a vector order parameter instead of 
0. scalar. Here we report a study of the equilibrium 
Qtate of spinor condensates in an optical trap. The 
freedom of spin orientation leads to the formation 
—of spin domains in an external magnetic field. The 
Structure of these domains are illustrated in spin do- 
C^ain diagrams. Combinations of both miscible and 
j^Jmmiscible spin components were realized. 

I A variety of new phenomena is predicted || [| |7| for 
fHspinor condensates, such as spin textures, propagation of 
O^pin waves and coupling between superfluid flow and atomic 
OSpin. To date such effects could only be studied in superfluid 
~+^He, which can be described by Bose-Einstein condensa- 
tion of Cooper pairs of quasi particles having both spin and 
terbital angular momentum gj. Compared to the strongly 
r^iteracting 3 He, the properties of weakly interacting Bose- 
Einstein condensates of alkali gases can be calculated by 
Ghean field theories in a much more straightforward and siru- 
JlJlfe way. 

^ i Other systems which go beyond the description with a 
Jingle scalar order parameter are condensates of two dif- 
arent hyperfine states of 87 Rb confined in magnetic traps. 
O&ecent experimental studies have explored the spatial sep- 
aration of the two components |l^| and their relative 
phas e [pTl| - Several theoretical papers describe their struc- 
ture | IX 11, 13, [l5| [Tg| [l7[ [l8| and their collective excitations 

Compared to these two-component condensates, spinor 
condensates have several new features including the vector 
character of the order parameter and the changed role of spin 
relaxation collisions which allow for population exchange 
among hyperfine states without trap loss. In contrast, for 
87 Rb experiments trap loss due to spin relaxation severely 
limits the lifetime. 

We consider an F=l spinor condensate subject to spin 



relaxation, in which two mp = atoms can collide and pro- 
duce an tof = +1 and an mp — — 1 atom and vice versa. 
We investigate the distribution of hyperfine states and the 
spatial distribution in equilibrium assuming conservation of 
the total spin. 

The ground state spinor wave function is found by min- 
imizing the free energy || 
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where kinetic energy terms are neglected in the Thomas- 
Fermi approximation which is valid as long as the dimen- 
sion of spin domains (typically 50 \xm) is larger then the 
penetration depth jll| (typically 1 \im). V is the trapping 
potential, n is the density, F is the angular momentum per 
atom, and E ze is the Zeeman energy in an external magnetic 
field. The Lagrange multiplier po accounts for the total spin 
conservation. The mean field energy in Eqn. (|l|) consists 
of a spin independent part proportional to c and a spin 
dependent part proportional to 02(F) 2 . The coefficients Co 
and C2 are related to the scattering lengths ao and a 2 for 
two colliding atoms with total angular momentum F to t = 
or F^t = 2 by c = Air h 2 a/M and c 2 = 4irh 2 Aa/M with 
a = (2a 2 + do)/3, Aa = (a 2 — ao)/3, and M for the atomic 
mass [^). The spin dependent interaction originates from 
the term c 2 i<i • F2 in the interaction of two atoms, which is 
ferromagnetic for c 2 < and anti-ferromagnetic for c 2 > 0. 

In the Bogoliubov approach the many-body ground state 
wave function is represented by the spinor wave function 



* if) = VW) C(r) = VW) (C+ (f) , Co (r) , C- (r)) , (2) 



where £+, Co, C- denote the amplitudes for the mp ■■ 
states, respectively, and |C| 2 = 1. 

The Zeeman energy E ze is given by 
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E ze = -E+IC+ 
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E_\C-\ 2 = E -p(F z )+q(F 2 ). 

(3) 

E + , Eq, E- are the Zeeman energies of the mp = +1,0, —1 
states, 2q = E + + £L — 2Eq is the Zeeman energy difference 
in a spinflip collision, and 2p = E- — E + . The Eq term can 
be included in the trapping potential V. The parameter p 
can be combined with the Lagrange multiplier pq to give 
p = p + p . 

In the following we determine the spinor which minimizes 
the spin-dependent part K s of the free energy: 



c{F) 2 -p(F z ) + q(F 2 ), 



(4) 



where c = c 2 n/2. The minimization of Eqn. (Q) for different 
values of the parameters c, p, and q is straightforward, and 
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is shown graphically in the form of spin-domain diagrams in 
Fig. 1. 

Experimentally, the values of c, p, and q can be varied 
arbitrarily, representing any region of the spin domain dia- 
gram. The magnitude (but not the sign) of the coefficient 
c is varied by changing the density n, either by changing 
the trapping potential, or by studying condensates with dif- 
ferent numbers of atoms. In this study, the axial length 
of the trapped condensate is more than 60 times larger 
than its radial size, and thus we consider the system one- 
dimensional, and integrate over the radial coordinates, ob- 
taining n — 2no/3 where no is the density at the radial 
center. This integration assumes a parabolic density profile 
within the Thomas-Fermi approximation. The value of q 
can be changed by applying a weak external bias field -Bo; 
q then corresponds to the quadratic Zeeman shift q = <jBq. 
The coefficient p arises both from the linear Zeeman shift 
and from the Lagrange multiplier po which is determined by 
the total spin of the system. For a system with zero total 
spin in a homogenous bias field Bq, pq cancels the linear 
Zeeman shift due to B , yielding p — 0. Positive (negative) 
values of p are achieved for condensates with a positive (neg- 
ative) overall spin. Finally, the coefficients can be made to 
vary spatially across the condensate. In particular, applying 
a field gradient B' along the axis of the trapped condensate 
causes p to vary along the condensate length. For a con- 
densate with zero total spin, p = pB' z where z is the axial 
coordinate with z = at the center of the condensate. Thus, 
the condensate samples a vertical line in the spin domain di- 
agrams of Fig. 1. The center of this line lies at p = 0, and 
its length is given by the condensate length scaled by fj,B' . 

The experimental study of spinor condensates required 
techniques to selectively prepare and probe condensates in 
arbitrary hyperfine states. Spinor condensates were pre- 
pared in several steps. Laser cooling and evaporative cooling 
were used to produce sodium condensates in the mp = — 1 
state in a cloverleaf magnetic trap p^| . The condensates 
were then transfered into an optical dipole trap consisting 
of a single focused infrared laser beam Q . Arbitrary popu- 
lations of the three hyperfine states were prepared using rf 
transitions Q|. After the spin preparation, a bias field Bq 
and a field gradient B' were applied for a variable amount 
of time (as long as 30 s), during which the atoms relaxed 
towards their equilibrium distribution, as shown in fig. 2. 

The profiles in Fig. 3 were obtained from vertical cuts 
through absorption images. They provide clear evidence of 
anti-ferromagnetic interaction. The spin structure is consis- 
tent with the corresponding spin domain diagram in fig. 1 a. 
Overlapping mp = ±l clouds as observed are incompatible 
with the assumption of ferromagnetic interaction. 

The strength c = (50 ± 20) Hz of the anti-ferromagnetic 
interaction was estimated by determining zi,, the location 
of the mp = to the mp — ±1 boundary, and by plotting 
p = pB' zi, versus the quadratic Zeeman shift q = qB^ as 
shown in Fig. 4. With n = (2.9 ± 0.5) x 10 14 cm" 3 the 
difference between the scattering lengths can be determined 
to a 2 -a = 3Aa = (3.5 ± 1.5)a B = (0.19 ± 0.08) nm where 
ob denotes the Bohr radius. This result is in rough agree- 



ment with a theoretical calculation of a 2 — clq — (5.5±0.5)a_e 
. The anti-ferromagnetic interaction energy corresponds 
to 2.5 nK in our condensates. Still, the magnetostatic (fer- 
romagnetic) interaction between the atomic magnetic mo- 
ments is about ten times weaker. It is interesting to note 
that the optically trapped samples in which the domains 
were observed were at a temperature of the order of 100 
nK, far larger than the anti-ferromagnetic energy. The for- 
mation of spin domains occurs only in a Bose-Einstein con- 
densate. 

Fig. 3 c shows a profile of the density distribution for a 
cloud at Bq = 20 mG and almost canceled gradient (B' < 
2mG/cm). No mp = region can be identified. The cloud 
was prepared with a small total angular momentum. Due to 
the almost-zero gradient and the non-zero angular momen- 
tum the cloud corresponds to a point in the shaded region 
in Fig. 1 a, rather than a vertical line with no offset as dis- 
cussed before with finite gradients and zero angular momen- 
tum. The different widths of the profiles are probably caused 
by residual field inhomogenitics. Fig. 3 c demonstrates the 
complete miscibility of the mp = ±1 components. 

For a homogenous two-component system the criterion 
for miscibility (immiscibility) is 

dab < {>)^Ja a ab JlJ, [if], [l8), when the mean field energy is 
parametrized as (2ttTi 2 /M)(n^a a + n 2 ab + 2n a riba a b)- Here, 
n a fi and a a ^ are densities and scattering lengths for the 
components a and b, and the scattering length a a t charac- 
terizes the interactions between particles a and b. In our 
spinor condensate with mixtures of the mp ± 1 components, 
we have a_i = a + i = a + Aa and a_i + i = a — Aa. Thus 
Aa > 0, like experimentally observed, implies miscibility. 
For a mixture of the mp = 1 and mp — components, we 
find do = a, a + i = a + Aa and ao+i = a + Aa, corre- 
sponding to immiscibility. For the 87 Rb experiments ||, [Io| 
it is not clear whether the two components are miscible or 
overlap only in a surface region due to kinetic energy pq] . 

In conclusion, Bose-Einstein condensates of sodium oc- 
cupying all three hyperfine states of the F = 1 ground state 
multiplet were optically trapped in low magnetic fields. The 
hyperfine states are coupled by spin-exchange processes, re- 
sulting in the formation of spin domains. We developed 
spin domain diagrams for both the anti-ferromagnetic and 
the ferromagnetic case, and showed that sodium has anti- 
ferromagnetic interactions, whereas the opposite case is pre- 
dicted for the 87 Rb F = 1 spin multiplet P4j| - All regions in 
the spin domain diagrams are accessible with our experimen- 
tal technique and thus any combination of the three hyper- 
fine components can be realized by applying small external 
magnetic fields. Of special interest for future work is the zero 
magnetic field case, where the rotational symmetry should 
be spontaneously broken. We observed both miscibility and 
immiscibility of hyperfine components. Thus the dynamics 
and possible metastable configurations j7| of two interpen- 
etrating, miscible superfluid components (mp = ±1) with 
arbitrary admixtures of an immiscible component (to p = 0) 
can now be studied. 
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Figure caption 1: Spin domain diagrams for spin-one 
condensates. The structure of the ground state spinor is 
shown as a function of the linear (~ p) and quadratic (~ q) 
Zeeman energies. Hyperfme components are mixed inside 
the shaded regions. Solid lines indicate a discontinuous 
change of state populations whereas dashed lines indicate 
a gradual change. The behaviour for q < is also shown al- 
though it is not relevant for this experiment. For c = 0, the 
Zeeman energy causes the cloud to separate into three do- 
mains with TTi p = +1,0, —1 and with boundaries at \p\ = q, 
as shown in b. For C2 ^ 0, the mean field energy shifts 
the boundary region between domains and leads to regions 
of overlapping spin components. In the anti-ferromagnetic 
case (a), the mp — component and the mp = ±1 com- 
ponents are immiscible (including the kinetic energy terms 
in Eqn. (Q) would lead to a thin boundary layer) and the 
boundary occurs at \p\ = q+c. For small bias fields, with q < 
c and \p\ < 2c, the mp = domain is bordered by domains 
in which mp =±1 components are mixed. The ratio of the 
mp = ±l populations in these regions does not depend on q, 
but is given by |C+| 2 /|C-| 2 = (2c+p)/(2c — p). In this region 
of small fields, the boundary to the mp — component lies 
at \p\ — 2^/cq. In the ferromagnetic case (c) all three com- 
ponents are generally miscible, and have no sharp bound- 
aries. Pure mp = domains occur for \p\ < J q(q — 4|c|) 
and pure mp = ±1 domains for |p| > q. Here, in contrast 
to the anti-ferromagnetic case, a pure mj? = condensate 
is skirted by regions where it is mixed predominantly with 
cither the mp = — 1 or the mp = +1 component. The con- 
tribution of the third component is very small (< 2%). In 
all mixed regions the mp = component is never the least 
populated of the three spin components. This qualitative 
feature can be used to rule out that F = 1 sodium atoms 
have ferromagnetic interactions. 

Figure caption 2: Formation of ground state spin do- 
mains. Absorption images of ballistically expanding spinor 
condensates show both the spatial and hyperfme distribu- 
tions. Arbitrary populations of the three hyperfme states 
were prepared using rf transitions (Landau-Zener sweeps) 
Q . At a bias field of about 40 G the transitions from mp = 
— 1 to mp = and from mp = to mp = +1 differ in fre- 
quency by about 0.9 MHz due to the quadratic Zeeman shift 
and they could be driven separately. The images of clouds 
with various dwell times in the trap show the evolution to 
the same equilibrium for condensates prepared in either a 
pure mp — state (upper row) or in equally populated 
mp = ±1 states (lower row). Between 5 s to 15 s dwell time, 
the distribution did not significantly change, although the 
density decreased due to three-body recombination. The 
bias field during the dwell time was Bq = 20 mG and the 
field gradient was B' = 11 mG/cm. These images were taken 
after the optical trap was suddenly switched off and the 
atoms were allowed to expand. Due to the large aspect ra- 
tio (typically 60), the expansion was almost purely in the 
radial directions. All the mean-field energy was released af- 
ter less than 1 ms, after which the atoms expanded as free 
particles. Thus, a magnetic field gradient, which was applied 
after 5 ms time-of-flight to yield a Stern-Gerlach separation 



of the cloud, merely translated the three spin components 
without affecting their shapes. In this manner, the single 
time-of-flight images provided both a spatial and spin-state 
description of the trapped cloud. Indeed, the shapes of the 
three clouds fit together to form a smooth total density dis- 
tribution. After a total time-of-flight of 25 ms the atoms 
were optically pumped into the F = 2 hyperfme state and 
observed using the mp — +2 to mp = +3 cycling transi- 
tion. This technique assured the same transition strength 
for atoms originating from different spin states. The size of 
the field of view for a single spinor condensate is 1.7 mm x 
2.7 mm. 

Figure caption 3: Miscible and immiscible spin do- 
mains. Axial column density profiles of spinor Bose-Einstein 
condensates are shown, obtained from time-of-flight absorp- 
tion images as in Fig. 2. The profiles of the mp = ±1 
components were shifted to undo the Stern-Gerlach sepa- 
ration. At low bias fields (Fig. 3 a), the mp = compo- 
nent was skirted on both sides by mp = ±1 components 
with significant mp = =pl admixtures thus demonstrating 
the anti-ferromagnetic interaction (also visible in Fig. 2). 
At higher fields (Fig. 3 b), the mp = ±1 components 
are pushed apart further by a larger mp = component 
and the mp = +1 admixtures vanish. They could not be 
resolved for quadratic Zeeman energies q > 20 Hz. The 
anti-ferromagnetic interaction leads to immiscibility of the 
mp =0 and the mp = ±1 components. The kinetic energy 
in this boundary region, which is small compared to the 
total mean-field energy, is released in the axial direction. 
Due to this axial expansion of the cloud in the time-of- 
flight and due to imperfections in the imaging system in- 
cluding the limited pixel resolution, the mp = to mj? = il 
boundary is not sharp. Fig. 3 c demonstrates the complete 
miscibility of the mp = ±1 components. The magnetic 
field parameters were Bq = 20 mG, B' = 11 mG/cm in a, 
B Q = 100 mG, B' = 11 mG/cm in b, and Bq = 20 mG, B' < 
2 mG/cm in c. 

Figure caption 4: Estimate of the anti-ferromagnetic 
interaction energy c. Plotted is the linear Zeeman energy 
p = IptB'zbl at the boundary between the mp — and 
mp =±1 regions versus the quadratic Zeeman shift q = <jBq. 
q = g 2 n B /l6h 2 v hfs = 278Hz/G 2 , and n = g s p, B /4h = 
700kHz/G, g s denotes the electron g-factor, v^fs the hy- 
perfme splitting frequency, and fi b the Bohr magneton. The 
solid line is a fit of the function \p\ = 2^/qc for q < c and 
\p\ = q + c for q > c. Extrapolating the linear part to zero 
bias field (dashed line) yields c — (50 ± 20) Hz. The data 
points at a given bias field represent p = fiB'zb for differ- 
ent gradient fields B' and thus mp = regions of different 
size. The scatter of these points is mainly due to a residual 
magnetic field inhomogenities, resulting in small deviations 
of the local gradient B' . The error bar represents the rel- 
ative error of all data points of 30 % in p and 5 % in q 
as estimated from the uncertainties in the magnetic field 
calibration. Furthermore the limited pixel resolution and 
contributions of the kinetic energy in the condensate to the 
axial expansion enhance the errors for the determination of 
Zb of small mp = regions. 
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